Dielectric constant of ice Ih and ice V: a computer simulation study.
نویسندگان
چکیده
In this work, we exploit an efficient algorithm for the sampling of hydrogen-bond networks in order to study the order/disorder transition and the dielectric properties of ice Ih. Our results show that all of the rigid point charge models studied, TIP4P/Ice, TIP4P/2005, TIP4P, and SPC/E, yield a very low dielectric constant as compared with experiment. The analysis of the polarization factor, G, reveals that the structure of ice Ih as predicted by the TIP4P family of models studied is very similar. For all such models, G approximately 2.5, the differences in the dielectric constant being given mainly by the different molecular dipole moments. The angular correlations in the SPC/E model are very different, however, yielding G approximately 1.8. Our study suggests that the dielectric anisotropy of ice Ih is very small and at any rate smaller than the statistical uncertainty of our results. We show that for sufficiently anisotropic structures, such as ice V, the method employed does indeed reveal a strong anisotropy. At 5000 bar and 180 K, we find that the principal dielectric constants of ice V according to the TIP4P/Ice model are 34, 164, and 76. The order/disorder transition could not be well characterized, because of very strong hysteresis. Whereas low temperature ordered phases were readily disordered at temperatures below 50 K, the high temperature phase remained disordered upon cooling. Nevertheless, our plots of dielectric constant versus temperature reveal a clear lambda-like shape and obey a Curie-Weiss law. All of the models studied show a higher stability of the antiferroelectric Pna2(1) phase than the experimental ferroelectric Cmc2(1) phase.
منابع مشابه
Dielectric constant of ices and water: a lesson about water interactions.
In this paper, the dielectric constant has been evaluated for ices Ih, III, V, VI, and VII for several water models using two different methodologies. Using Monte Carlo simulations, with special moves to sample proton-disordered configurations, the dielectric constant has been rigorously evaluated. We also used an approximate route in which proton-disordered configurations satisfying the Bernal...
متن کاملDielectric constant and proton order and disorder in ice Ih: Monte Carlo computer simulations
متن کامل
Simulations of proton order and disorder in ice Ih.
Computer simulations of ice Ih with different proton orientations are presented. Simulations of proton disordered ice are carried out using a Monte Carlo method which samples over proton degree of freedom, allowing for the calculation of the dielectric constant and for the examination of the degree of proton disorder. Simulations are also presented for two proton ordered structures of ice Ih, t...
متن کاملProton ordering dynamics of H2O ice.
From high precision measurements of the complex dielectric constant of especially prepared samples of H2O, we identify the onset temperatures of the phase transition into and out of ice XI from ice Ih to occur at TIh-XI = 58.9 K and TXI-Ih = 73.4 K. For D2O, TIh-XI = 63.7 K and TXI-Ih = 78.2 K. A triple point is identified to exist at 0.07 GPa and 73.4 K for H2O and 0.08 GPa and 78.2 K for D2O ...
متن کاملELASTIC CONSTANTS OF ICE Ih, UP TO 2.8 KBAR, BY BRILLOUIN SPECTROSCOPY
The technique of high resolution Brillouin spectroscopy has been used to determine accurately the hydrostatic pressure dependence of the elastic constants of ice Ih in the full pressure range of phase stability 0 2.8 kbar. The percentage changes in the adiabatic elastic constants c,,, cI2. cI3, cS3, c,, and bulk modulus, over the pressure range, are 3.0, 7.2, 8.4, 2.8, -1.5 and 5.6 percent per ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- The journal of physical chemistry. B
دوره 114 18 شماره
صفحات -
تاریخ انتشار 2010